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We study the binding of a holon to a doublon in a half-filled Hubbard model as the mechanism of the zero- 
temperature metal-insulator transition. In a spin polarized system a single holon-doublon (HD) pair exhibits a 
binding transition on a 3D lattice, or a sharp crossover on a 2D lattice, corresponding well to the standard Mott 
transition in unpolarized systems. We extend the HD-pair study towards non-polarized systems by considering 
more general spin background and by treating the finite HD density within a BCS-type approximation. Both 
approaches lead to a discontinuous transition away from the fully polarized system and give density correlations 
consistent with numerical results on a triangular lattice. 

PACS numbers: 71.27,+a, 71.30.+h, 71.10.Fd 


I. INTRODUCTION 


The mechanism of the Mott metal-insulator transition 
(MIT) 1 , i.e., the transition in a half-filled band induced by 
the electron-electron repulsion, has been a major theoretical 
and experimental challenge for the last few decades 2 . The 
primary signature of the MIT (at zero temperature T = 0) 
is the appearance of the charge gap A c and vanishing of the 
charge susceptibility \c on the insulating side. A number of 
distinct mechanisms for the transition have been proposed. 
These include those due to Brinkman and Rice 3 (where the 
quasi-particle weight in the metallic phases approaches zero 
as the transition is approached), Hubbard (where vanishing 
of the charge gap occurs when the upper and lower Hubbard 
bands overlap), or Dynamical Mean-Field Theory (DMFT) 4 
which combines both these features. Powerful numerical ap¬ 
proaches have been applied to prototype models such as the 
Hubbard model, in particular methods emerging from DMFT 
and variational Monte Carlo approaches 5 . The physics of the 
MIT is complicated and challenging due to the possibility of 
coexisting antiferromagnetic (AFM) spin orderings within the 
Mott insulator, depending on the underlying lattice. The ef¬ 
fect of the competing spin order on the MIT is expected to 
be less pronounced in non-bipartite lattices, as in the Hub¬ 
bard model on a triangular lattice 5 ’ 9 , of relevance to organic 
charge-transfer salts 10 ’ 13 . Recently also fullerenes 14,15 have 
become of interest as an example of the MIT in 3D frustrated 
lattices. 

In this paper we show that within the half-filled single¬ 
band Hubbard model the MIT can be understood as a tran¬ 
sition from an unbound to a bound state of doubly occupied 
sites (doublons) and empty sites (holons). This view has been 
advocated long ago 1,6,16 ’ 18 , but made more explicit in more 
recent variational 5,19 , field theory 20 , and slave-boson 21 ap¬ 
proaches. The concept sketched for a triangular lattice in the 
right panel of Fig. 1 is simple: in a half-filled band one can 
consider the holon-doublon (HD) pairs as basic charge excita¬ 
tions relative to the reference state with singly occupied sites. 
Within the Mott insulator at large Coulomb repulsion U > U c 
the HD pairs (at T = 0) are bound while in a metal at U < U c 
pairs are not bound and holons and doublons separately con¬ 


tribute to electrical current and to the finite charge stiffness 
D c > 0. Such an approach to the MIT becomes even more 
reasonable if one considers a polarized half-filled band, i.e. 
the system with nearly saturated spin magnetization m 1/2 
due to, e.g., the presence of a uniform magnetic field. In this 
case one deals with a simplified problem of vanishing (but 
non-zero) concentration of HD pairs 22 , which can be treated 
exactly on any lattice model. At half-filling introducing the 
magnetization (|m| < 1/2) as an additional degree of free¬ 
dom opens the possibility to explore the ground-state phase 
diagram of the Hubbard model for any particular lattice. 

In particular, this approach allows for a novel path to the 
challenging regime of m ~ 0, as well as a benchmark and an 
alternative interpretation of the MIT. In Fig. 1 we anticipate 
the result for the Hubbard model on a triangular lattice, which 
will be our focus later on. Extending the study of the MIT to 
Hubbard systems in large magnetic fields 23,24 is of experimen¬ 
tal relevance, e.g., for organic conductors in a magnetic field 25 
and for cold fermions in optical lattices. We should mention 
here the relation to several studies of the attractive Hubbard 
model 26 " 29 , which can be mapped to the repulsive model in a 
magnetic field 30 , provided that the lattice is bipartite. 
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Figure 1. (Color online) Phase diagram of the Hubbard model on a 
triangular lattice at half filling as calculated from the BCS-type ap¬ 
proximation and exact diagonalization (Lanczos) results. The Lanc- 
zos results are for magnetization m / 0 calculated at Nd holon- 
doublon pairs and system sizes N with values Nj/N = 2/100, 2/36 
and 3/36, while result for m = 0 is from Ref. 8. Right panels show a 
schematic unbound (bound) holon-doublon pair in the metallic (insu¬ 
lating) phase on a triangular lattice with a non-zero spin polarization. 
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The paper is organized as follows: In Sec. II we employ an 
analysis of the single HD pair in a spin polarized background 
m ~ 1/2 . Concentrating on non-bipartite lattices, in particu¬ 
lar the triangular lattice, this simplified case gives the critical 
(or crossover) U c for the MIT transition and also the depen¬ 
dence of the charge gap A c ([7), and double occupancy D that 
are quite close to the numerical results for the unpolarized sys¬ 
tem. The analysis is extended towards the non-polarized sys¬ 
tem in two ways. On one hand we consider in Sec. Ill the 
motion of a HD pair within a more general spin background, 
for which the MIT is a discontinuous one. In an alternative 
approach we deal in Sec. IV directly with a finite density of 
HD pairs within a BCS-type approximation 6 - 26,31 , giving a fair 
agreement with numerical finite-system results. Conclusions 
and implications are presented in Sec. V. 


II. SINGLE HOLON-DOUBLON PAIR BINDING 


We consider the single-band Hubbard model with nearest- 
neighbor (n.n.) hopping on a general lattice, 

H = -tJ2 ( c ]a c i* + c\ a c ja ) + UJ2 n^nn- (1) 

(ij)a i 

We restrict our study to a half-filled band and T = 0. To 
analyze the MIT we consider the system with a finite magne¬ 
tization m = St ot /N, where Sf ot = JV S'? is the total spin, 
and N the number of lattice sites. Clearly, the case with max¬ 
imum m = 1/2 is a ferromagnetic (FM) band insulator, but 
the problem is nontrivial for any \m\ < 1/2. 

We first consider the problem within the sector S* ot = 
N/2 — 1, i.e. a single HD pair, (we note that a similar case 
has been considered also by Kohn 16 ). The problem within the 
model (1) is that of two particles with total momentum q, 


\r q ) = N-V 2 Y,j ql Wl,l+r), Wlm ) = clc mt \FM), 

l 

( 2 ) 

and the onsite attraction —U, leading to the gap equation for 
the pair energy E q 22 , 


1 = - 


U 

N 


E 


i 

E q -U -r) k ’ 


Tjk — €q—k - (3) 


Here e/- is the single-electron band energy. Eq. (3) yields the 
energy E q of the bound HD pair at given q provided that E q < 
U + rnirifc [r/fc]. In the opposite case, the pair is unbound and 
states form a continuum. For a general lattice the minimum 
for the bound HD pair is Eq = min[f? r / = E qo whereby qo is 
straightforwardly determined for standard cases. 


A. Hypercubic lattices 

In a D-dimensional hypercubic lattice one gets, due to the 
particle-hole symmetry, qo = 7r(l,...,l) and t qo -k = — fifc 
so that ? 7 fc = —2efc. The continuum of states has the lower 


edge at E\ = U — 2 zt, z being the number of n.n., so that the 
charge gap 


A c — Ei — Eq — U — 2 zt — Eq. (4) 

In the insulating regime ( U t) it follows that Eq ~ 

—4zt 2 /U, and so A c ~ U — 2 zt 1 . Finite U c ~ 2zt > 0 ap¬ 
pears in a 3D case. On the other hand, in 2D there is strictly no 
transition due to the singular (step-like) single-electron den¬ 
sity of states (DOS) V(e) at the band edge. The latter in con¬ 
nection with Eq. (3) generally leads for 2D to U c = 0. Still 
there is a sharp crossover in A C (U ~ U*) with A c being ex¬ 
ponentially small for U < U*. The same effect is present 
also within the triangular lattice discussed further-on with the 
result for A c presented in Fig. 2. 

The m ~ 1/2 case should be distinguished from m = 0, 
where, e.g., on a square lattice U c = 0 due to Fermi surface 
nesting 32 - 33 . Here, we should mention a related limit for the 
z —» oo Bethe lattice, relevant for DMFT studies, where the 
central input is the DOS 

D(e) = (2/-kW)VW 2 -e 2 , (5) 

with W being the effective half-bandwidth. In this case one 
finds analytically that U c = W, serving as a nontrivial bench¬ 
mark for numerical studies 29,34 . 



U/r 

Figure 2. (Color online) Charge gap A c vs. U for the Hubbard 
model on triangular lattice. The results shown are calculated for a 
single holon-doublon (HD) pair at m ~ 1/2 (thick blue line), a sin¬ 
gle HD pair in an unpolarized spin background m = 0 within the 
retraceable-path approximation (thin black line), and numerically for 
m = 0 (red line with points). 


B. Triangular lattice 

Let us now turn to the non-bipartite lattices, where a widely 
studied case is the Hubbard model on a triangular lattice with 
n.n. hopping. E q reaches a minimum value at qo = (47t/3, 0). 
Since we are dealing with 2D DOS there is strictly no MIT but 
rather a sharp crossover at U* ~ 6 1, well visible in A C (U), 
presented in Fig. 2. We can judge the relevance of this exact 
result for to ~ 1/2 by comparing it with numerical results 
for the case m = 0, obtained by using the finite-temperature 
Lanczos method (FTLM) on N = 16 sites (taken from Fig. 2 
in Ref. 8 and multiplied by 2 due to different definitions of 








3 


the gap). In spite of quite different spin background. Fig. 2 
shows a fair correspondence between the single-HD result (at 
least for U/t > 10) and A C (U) at to = 0. Clearly, one should 
take into account that at m ~ 1/2 the gap A C (U) exhibits 
only a crossover at U* (with exponentially small gap for U < 
U*) while strictly U c = 0 for m = 1/2. One should bear 
in mind that also other numerical results for triangular lattice 
show considerable uncertainties, see e.g. the compilation in 
the Supplementary material for Ref. [8], 

For a single HD pair it is straightforward to evaluate also 
the actual doublon density D, defined by the probability that 
the holon and doublon are not on the same site, leading to the 
explicit expression. 


D = l- 


(~y _-_ 

\N ^ (E 0 -U-r ) k ) 2 


-1 


( 6 ) 


We calculate double occupancy as D = D /2 to have correct 
normalization per site in order to compare with D standard for 
m = 0. We then compare it with other results for m < 1/2 in 
Fig. 3. 



Figure 3. (Color online) Double occupancy D (per double number 
of pairs 2 Nd) vs. U. The results shown are calculated for m ~ 1/2 
with a single HD pair (denoted 1HD), for a single HD pair within 
the RPa (denoted RP approx.), with FTLM and taken from Ref. 8 
(denoted FTLM), with exact diagonalization at T = 0 and various 
m (denoted ED), and with BCS approximation for particular m (de¬ 
noted BCS). 


III. HOLON-DOUBLON PAIR IN A GENERAL SPIN 
BACKGROUND 

An extension from m ~ 1/2 towards to = 0 can be made 
by keeping a single HD pair, but formed in a more general 
spin background. The reference state should be first found 
within the spin Heisenberg model describing the Mott insulat¬ 
ing phase. Being primarily interested in frustrated lattices we 
assume that spin correlations are of short range. In such a case 
the retraceable path approximation (RPa) is appropriate 35 , ex¬ 
cluding the HD motion along loops. Within the RPa one de¬ 
scribes a HD pair created/annihilated on the same site with the 


local Green’s function, 

G(u>) = -——E(w) = 2azt 2 G2{io)i (7) 

UJ — 2-i yiO) 

where the recombination hopping is renormalized by the spin 
correlation factor 

a = <^ol(^ — 2S, • Sj)|$ 0 ) = ^ — 2/x. (8) 

Here i. j denote nearest neighbors. Further hoppings of either 
holon and doublon are taken into account with G 2 (w). As¬ 
suming that both holon and doublon hop independently along 
the retraceable paths and neglecting spin string energy, G 2 {uj) 


, ( 9 ) 

J J UJ -U -ex -62 

where the DOS within the RPa is given for the correspond¬ 
ing Bethe lattice, Eq. 5, with W = 2 y/z — It. 36 . The bound 
solution {Eq = 'E(Eq)) is stable provided that it is below 
the continuum of unbound HD-pair states at E\ = U — 2W 
derivable from Eq. (9). The corresponding energy difference 
A c = Ex — Eq is interpreted as the charge gap A c rele¬ 
vant for the MIT. In contrast to the single HD pair in the FM 
background, Eq. (3), the MIT is a discontinuous one due to 
energy level crossing of bound and unbound HD-pair solu¬ 
tions. It should be mentioned that direct comparison with the 
m ~ 1/2 is not feasible since loop trajectories would be re¬ 
quired to represent fully the latter case. In Fig. 2 we plot the 
resulting A C (U) taking /i = —0.182 for a as appropriate for 
the Heisenberg model that holds for U t™°. It is ev¬ 
ident that the RPa result is closer to the numerical one for 
to = 0, which is obtained by FTLM and shows small finite 
size effects 8 . This agreement gives support to the RPa for tri¬ 
angular lattice. Within the same method it is straightforward 
to evaluate also the effective double occupancy (per HD-pair) 
D evaluated from Eq. (7) via the weight of bound-state pole 

( 10 ) 

or from dE q/OU. We present the results in Fig. 3, where it 
can be seen to compare favourably with the FTLM result for 
to = 0 . 


IV. FINITE HOLON-DOUBLON PAIR DENSITY 

An alternative approach to the case of a finite HD pair den¬ 
sity Nd/N with maximal number of doublons Nd and holons 
Nh = Nd, and m = 1/2 — Nd/N, is to start from a reference 
FM state (all spins up) and to define doublon d\ = cj, and 

holon h\ = Cj-f- operators, and rewrite the repulsive Hubbard 
model, Eq. (1), as an HD attractive one 6 - 30,41 , 

H = + e h k h\h k ) + UN d -UY n diU hi . (11) 

k i 
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Here, ef = —e k = ek- This should be distinguished from 
the standard negative-!/ model with fermions having the same 
dispersion. Due to attraction between holons and doublons it 
is natural to approximate the ground state with a mean-field 
BCS-type approach, since it captures HD pair binding and 
their coherent superposition. 

A. BCS-type approach 

As in several previous works 6,26 ’ 31,42 we approximate the 
ground state of the collective HD pair state with a BCS- 
type wavefunction, representing the coherent superposition of 
bound HD pairs, 

|'T°) = + v k b\)\FM) (12) 

k 

where b\ = d\h/ q _ k is the HD pair creation for given q and 

C ulvl) = (l±( k /E k )/2 . (13) 

Here, \FM) is the reference FM state of all spins up as already 
used in Eq. (2). </k = r/fe/2 — p (with rjk from Eq. (3)) is 
the HD pair dispersion, p = {p d + ph )/2 is the effective 
pair chemical potential, Ek = \/Cf + A 1 2 the quasi-particle 
excitation energy and A the pairing order parameter. The BCS 
ground state energy is then given by 

Ebcs = »'*• “ + Jj l A ! 2 + W 1 - *d) + 2 P\ N d> ’ 

k 

(14) 

where we used the abbreviation n d = Nd/N. We note that 
j'Fo), Uk, Vk, Ek, Ck and A depend also on q, and we omit this 
from our notation for clarity. At fixed rid (or m = 1/2 — n d ) 
one still has to solve self-consistently equations for p and A, 

1 v- 2 1 T 1 ... 

"**»?"*■ 1 = (5) 

In addition, q that gives the lowest energy Ebcs has to be 
determined. Nevertheless, we consider here only q = qo = 
(47t/3,0), which is clearly the solution for m ~ 1/2 and also 
away from m = 0 6 . 

A does not have the meaning of the charge gap 26 like in the 
usual weak-coupling BCS theory for superconductors. For 
the MIT one needs to compare Ebcs at a given fid with the 
metallic state or HD liquid. The energy of the metallic state 
Em is obtained from Eq. (14) by putting A = 0 and filling 
the doublon and holes states up to pd and ph, respectively. 

The MIT appears as a crossing of Ebcs and Em at U = U c 
or the appearance of positive condensation energy E c = 
Em — Ebcs with increasing U. The transition is of first or¬ 
der due to jump of A at this crossing. Within BCS one can 
calculate also the double occupancy D (normalized to 2N d ) 
as 


It is important to realize that within a BCS approximation 
one can reproduce the exact single HD pair results. We note 
that within the dilute limit, fid <C 1, where also (k 3> A, 
Eqs. (15) reduce to the form for a single HD pair, Eq. (3), 
on any lattice. The relation is derived by identifying p with 
the single HD pair energy E = U + 2p obtained within the 
rid <C 1 limit of Eq. (14). In this case, it also follows 

Ck ~ U/2, A ~ 2C,k\/fid ~ Uy/Wd <C t, (17) 
in analogy to bipartite lattices 26,30 . 

B. Comparison with numerical results 

Within BCS one can also calculate the double occupancy 
D, Eq. (16), and we compare the result in Fig. 3 with the 
exact diagonalization (ED) result obtained with the Lanczos 
method 43 on a finite cluster for Nd/N = 1/12 (Nd = 3, 
N = 36) and find good agreement. ED results are shown 
also for Nd = 2 HD pairs on N = 36 and N = 100 sites, 
respectively. 

The validity of the BCS approximation at low HD densi¬ 
ties Nd/N 1/2 can be tested also via density correlation 
functions. Of particular interest are the doublon-holon and 
doublon-doublon (or holon-holon) density correlations 

C d h(l'') — ^ '^ (fldirihd+r) •> CddiX*) — ^ , 

i i 

(18) 

where n di = n di {l-n hi ), n hi = n hi {l-n di ) are ’true’ (pro¬ 
jected) doublon and holon operators, respectively. We present 
the comparison between results from ED at N d = 2 and the 
BCS function also projected onto N d = 2 HD pairs. Fig. 4 
demonstrates the effective HD attraction via C d h and the re¬ 
pulsion between doublons via C dd - The agreement between 
ED and the BCS approximation is remarkable within the in¬ 
sulator regime U/t = 6,8 > U c /t. Regarding correlations, 
the main signature of the insulating (HD binding) regime is 
that C d h and C dd converge (according to the BCS approxima¬ 
tion exponentially) at larger r. In contrast, in the metallic case 
the approach is expected to be much slower (algebraic) and 
indeed for U = 4f the C d h remains larger than C dd even for 
the largest r in a given system indicating that the holon is not 
bound to the doublon. 

Finally let us discuss the phase diagram of the Hubbard 
model on a triangular lattice, presented already in Fig. 1. In 
the limit m ~ 1/2 the BCS insulating state is stable for any 
U > 0. On the other hand, for any m < 1/2 the condensa¬ 
tion energy E C (U) is vanishing at finite U = U c . For com¬ 
parison we present in Fig. 1 also ED results, where the MIT 
transition is monitored by the discontinuity of the double oc¬ 
cupancy D(U). Value of U c (m = 0) is taken from FTLM 
result presented in Ref. 8. 

V. CONCLUSIONS AND IMPLICATIONS 


1 - n d _ A 2 

2 2 n d U 2 ' 


( 16 ) 


The aim of this work is to present a theory of the Mott MIT 
within a half-filled Hubbard band starting from a spin polar- 
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Figure 4. (Color online) Doublon-holon Cdh and doublon-doublon 
Cdd correlations vs. distance r for different U/t evaluated for two 
HD pairs, Nd = 2, on N = 64 sites. Results obtained by using the 
full ED (lines) can be compared to the ones evaluated within the BCS 
wavefunction for U/t = 6,8 (points). 


ized system. The advantage of such an approach is that one 
can systematically follow the MIT by reducing the polariza¬ 
tion by which the behavior evolves from the exactly solvable 
(but nontrivial) single HD-pair binding at m ~ 1/2 into a col¬ 
lective HD-pair phenomenon going towards m ~ 0. Our con¬ 
sideration of m < 1/2 requires some approximations (RPa 
and BCS) but we show that they produce meaningful and quite 
reliable results also for m ~ 0. Within these approximations 
the MIT at m < 1/2 is discontinuous, and the following sim¬ 
ple scenario is suggested by the BCS. When U is reduced from 
large values where the density and size of HD pairs are small, 
the pairs start to increase in size and the kinetic energy de¬ 
creases. When the pairs start to overlap and due to the hard 
core repulsion between them (see Cdd in Fig- 4) the kinetic 
energy cannot be further reduced, the metallic state becomes 
preferred, leading to a discontinuous MIT. 

Concentrating on a triangular lattice we have shown that the 
charge gap A C (U) for a single-HD pair is quite similar to that 
for m ~ 0, bearing in mind that A C (U) at m ~ 1/2 appears 
only as a crossover. Similar agreement is seen for the (prop¬ 
erly normalised) double occupancy D(U). When we consider 
the agreement with the numerical results, we should point out 
that a part of the uncertainties are due to the approximations 
(RPa and BCS) used. Moreover, numerical benchmarks for 
the MIT (available in the literature primarily for m ~ 0) in 
the triangular lattice are even more scattered (see the compi¬ 
lation in Ref. 8), so it is at this stage hard to better quantify 
our HD-pairing results. 

On the other hand, it is evident that presented RPa and BCS 
approaches have some restrictions when considering unpolar¬ 
ized m ~ 0 regime. Both approximations deal primarily with 
the bound pairs, i.e. the insulating phase, while the metallic 
side of the MIT is considered only on the simplest level, which 
can to our opinion underestimate the critical U c of the MIT at 
m ~ 0. Within RPa we did not take into account the influ¬ 
ence and competition of different spin orders close to the MIT, 
which should be less severe for frustrated lattices. With in¬ 
creasing concentration of HD pairs the interplay of close pairs 
is not taken into account within the RPa approach. However, 
it should be stressed that the criterion for the latter inter-pair 


effects is not directly Nd/N, but rather D which measures the 
density of ’real’ doublons and holons. In this respect, we no¬ 
tice that at ra ~ 0 D ~ 0.1 <C Nd/N is small even at the 
MIT, as evident for the triangular lattice, e.g. in Fig. 3. 

Concerning the BCS approximation it has two essential in¬ 
gredients, important to describe the holon-doublon (HD) insu¬ 
lator at finite HD density fid: a) BCS wavefunction represents 
the coherent superposition of bound HD pairs, b) it takes into 
account the Pauli exclusion principle at arbitrary fid . The lat¬ 
ter clearly leads to proper (repulsive - like) fermionic correla¬ 
tions Cdd and Chh- The correlations Chd are also taken on the 
two-body level. Still, it is evident that all correlation effects 
are not taken into account within the BCS wavefunction. The 
BCS approximation is also better justified for more extended 
pairs, while its non-straightforward limitations are similar to 
the strong-coupling BCS (see Ref. 26). BCS approximation 
also neglects quantum fluctuations of the gap order parameter, 
which could destabilise the insulating phase in some regime. 
It in addition yields the condensation of HD pairs, which is 
not necessary for MIT. It therefore remains to be checked to 
what extent are the BCS results, as e.g. Chd, quantitatively 
reproduced for fid 0. This is, however, a challenging task 
also for the present status of numerical methods. The valid¬ 
ity of the BCS phase diagram in Fig. 1 evidently depends on 
the quality of the BCS wavefunction, but also on the assumed 
simple form of the reference metallic state energy Em- Still, 
the sharp saturation of U c (m ~ 0.5) has (at least partly) the 
same origin to 2D singularity DOS, leading also to the very 
sharp crossover stemming from Eq.(3) and showing up in the 
A C (U) in Fig. 2. While further nearly constant U c {m) emerg¬ 
ing from BCS equations (calculated at fixed q = qo) has no 
evident explanation, it appears that only gross features of the 
band matter, as is also the case e.g. for A C (U > U* ) in Fig. 2. 

In our study we concentrated on the charge sector, being the 
driving mechanism of the MIT. It is however evident that away 
from m ~ 1/2 in the insulating phase the spin long-range or¬ 
der, e.g. the spiral (S q ) or an AFM order, can emerge. This 
appears to be the case also for the extensively studied frus¬ 
trated triangular lattice, at least for U ^ 7 - 37 - 38 - 44 - 51 . While 
it is evident that for bipartite lattices the nesting mechanism 
leads to the AFM order at m = 0 for any U > 0, the question 
of possible spin order in frustrated lattices remains challeng¬ 
ing, even for numerical methods. Still we should point out that 
the spin order is a secondary phenomenon for the MIT since 
it involves only exchange energies of the order of J = 4t 2 /U, 
whereas the charge sector leads to the MIT due to the HD at¬ 
traction U. We note that HD binding is also prerequisite for 
the description with Heisenberg type models, which are typi¬ 
cal starting points for the discussions of spin ordering. 

The character of the MIT transition at m ~ 1/2 depends 
on the lattice dimensionality. While in the 2D lattice the sin¬ 
gle HD-pair exhibits only a rather sharp crossover, the bind¬ 
ing transition becomes well defined and continuous in 3D lat¬ 
tices. The most intriguing is the f.c.c. lattice with U c ~ 12f 
being also of experimental relevance to fullerides 15 and a nat¬ 
ural extension of the presented work. Finally, we comment 
on the relation to other scenarios (or interpretations) of the 
Mott MIT. We are close to the original Mott proposal 1 for the 
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charge gap as A c = U — 2 zt. Similar closing of the gap be¬ 
tween Hubbard bands is incorporated also in the DMFT theo¬ 
ries. Still, the DMFT approach, as well as the Brinkman-Rice 
picture, is based on single-particle properties (the one-electron 
Green’s function) and the relation to the HD binding, a two- 
particle phenomenon is less evident. In any case, the relation 
of our HD-pair binding approach to existing theories of the 
MIT should be further explored. 
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